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Historically multi-item inventory control has been mod- 
eled by assuming that each item can be treated indeoendentlv 
in a variable cost minimization formulation. In this paper 
independence between items is not assumed. Constraints on 
total system operating characteristics create inter-item 
dependencies. Optimal policies are determined from a goal- 
constraint formulation. This is done without reliance upon 
unknown parameters such as order cost and carrying cost which 
the traditional theory leans on heavily. A group of models 
are presented, with necessary and sufficient conditions for 
optimal solutions provided for each. In addition, solution 
algorithms are indicated for the major models. An algorithm 
for verification of sufficiency conditions for a non-convex 
objective function is also provided. 
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CHAPTER I 



INTRODUCTION 

1 . 1 BACKGROUND 

The great majority of research in the field of inventory 
theory has been associated with the treatment of a single 
item. The object of such research is to determine decision 
rules which tell us how much to order and when to order to 
minimize variable costs. In order to implement these deci- 
sion rules one must know, or at least be able to estimate 
the order cost, holding cost, and shortage cost. Another im- 
plicit assumption is that investment capital is available to 
buy and hold the optimal order quantity and safety level. 

In the military the situation is somewhat different. 

The cost parameters mentioned previously are difficult, if not 
impossible, to estimate. Furthermore, they may not always be 
meaningful. An example of this is the inventory carrying cost. 
A major portion of the carrying cost is the opportunity cost 
related to alternative uses of the investment in stock. With- 
in the -Navy, the budget funds allocated to stock cannot be 
used for any other purpose. In addition to this problem, each 
military inventory is under severe budget restrictions. 

In order to understand these restrictions, some discussion 
of the budget process is necessary. The manager of a Naval 
Supply activity receives funds from two sources. The first 
source of funds is referred to as Operations and Maintenance 
funds. As the name suggests, these funds support the opera- 
tions of the activity which include salaries, providing for 
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maintenance of buildings, purchase of minor equiuments , arid re- 
plenishment of material necessary to the operation of the ac- 
tivity. These funds cannot be used for the acquisition of 
items for general stock. The second budget source provides 
funds for acquisition of items for general stock. This budget 
can be in terms of either a limit on average investment or a 
procurement limitation. Moreover, these resources cannot be 
transferred to the Operations and Maintenance budget. Thus 
one finds within the military existing constraints on funds 
available to purchase stock or those available to initiate 
and process orders and receipts. 

These differences between the military situation and the 
assumptions made in classical theory have suggested some of 
the questions to be explored in this study. 

These questions may be stated in a more precise form as 
follows : 

1. Is the single item theory sufficient to provide 
decision rules for the multi-item inventories in the 
Navy? 

2. Is the assumption that the order cost, holding 
cost, and shortage cost are known parameters reason- 
able or even necessary? 

3. Does the addition of budget constraints to the 
problem require a different basic model? 

It is proposed that a more adequate approach is that of a multi- 
item inventory model based upon an investment limit and an order 
constraint. 
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Depending upon the situation, the military inventory man- 
ager attempts either to minimize shortages or to minimize 
weighted shortages. In the first case the tidewater or first 
line stock point wants to minimize the number of times a cus- 
tomer is turned away due to a lack of material at that activ- 
ity. In the second instance, weighted shortages are minimized 
because some projects are more essential from a military 
standpoint than others. Thus, for the military, the cost of 
being out of an essential part is greater than the cost of be- 
ing out of an ordinary part. 

The above arguments have motivated the construction of a 
multi-item inventory model which attempts to satisfy one of 
the above goals subject to a constraint on average investment 
as well as a constraint on the number of orders per year. 

In order to have other than a trivial solution, every 
inventory model must have some built-in control to prevent 
either an infinite order quantity or, on the other extreme, 
an infinite number of orders. For instance, in the Economic 
Order Quantity inventory model the order cost prevents order- 
ing after each demand while the holding cost prevents an in- 
finite amount of inventory on hand. The optimal solution 
then achieves a trade-off between the order cost and the 
holding cost. 

Now we observe that a constraint on the average inventory 
would achieve the same purpose as the holding cost and a con- 
straint on number of orders would achieve the same purpose 
as the order cost. Moreover, the budget constraint and the 
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order constraint are more realistic since they are relatively 
easy to obtain and the inventory manager, who cannot be ex- 
pected to understand the mathematical theory intimately, has 
a better intuitive feeling for these constraints than he would 
have for the order cost and holding cost parameters. 



1.2 CURRENT NAVY PROGRAMS 

Currently the Navy has two programs which provide deci- 
sion rules for stock points and inventory control points. 

The first program which governs the management of retail 
material of Navy stock points is known as the Variable Opera- 
ting and Safety Level (VOSL) program. The VOSL program is a 
model of the following situation. The stock point receives 
an Operations and Maintenance allotment which, among other 
things, fixes the personnel ceiling at that activity. This 
in some sense limits the maximum number of orders which can 
be written. In addition, the stock point receives an invest- 
ment limit which cannot be exceeded. The stock point manager 
then strives to minimize some function of shortages while re- 
maining within the above constraints. In the case of the 
stock point we have essentially a single warehouse system. 

The second program which governs the purchase of Navy 
wholesale material is known as the Uniform Inventory Control 
Point program. However, in the case of the inventory control 
point, we have a multi-echelon, multi-warehouse system. On 
the ICP level, the inventory manager is given both a procure- 
ment and personnel budget. In this instance, the manager again 
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attempts to minimize shortages within the multi-echelon frame 
work. On the stock point level a similar situation exists 
with the exception that the inventory manager attempts to min 
imize some function of shortages with respect to his single 
warehouse system. 

In each program the manager attempts to minimize some 
function of shortages while subject to some constraint on in- 
vestment and a constraint on orders. This observation leads 
to the formulation of the models discussed in Chapter 2. 
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CHAPTER 2 



PROBLEM FORMULATION 

2.1 BASIC FORMULATION OF THE PROBLEM 

As suggested by the introduction, some other formulation 
of the inventory problem seems desirable. Such a formulation 
would not be based on minimizing variable costs in the classic 
sense, but rather on minimizing some more reasonable objec- 
tive. For this study the objective will usually be some func- 
tion of weighted units short. 

In general, the basic problem can be stated as 

min f(X,Y) 

subject to: 

gj(X,Y,A) £ 0 j = 1,1,. ..,m . 

where X and Y represent vectors of decision variables and A 
represents a vector of known parameters. Throughout this 
paper we will indicate vectors by non-subscripted capital let- 
ters. The objective function f(X,Y) can be expected units 
short, weighted units short or some other function of short- 
ages. The constraints can represent limitations on the ex- 
pected number of orders and restrictions on average invest- 
ment or a procurement budget. 

2.2 CONTINUOUS REVIEW FORMULATION 

If one is given an inventory of n items, an average in- 
vestment limit, and a limit on the number of orders then the 
following continuous review problem can be formulated. 
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2.2.1 Constraint Formulation 



Hadley and Whitin [1] show that under continuous review 
the average on hand quantity (m^) for the ith item is 



r . + — ^ - y . = m. , 
1 2 1 ' 



( 2 . 1 ) 



where 

r^ = the reorder point for the ith item, 

= the order quantity for the ith item, 

^i = /q X fi (x,T^) dx, 

f^(x,T^) is the density of lead time demand for the ith item, 
and r^ and (i = 1,2,* ‘'n) are the decision variables. 

To convert (2.1) to average on hand investment per item 



we multiply by the price of the ith item (c^) : 



Qi 

Ci(r^ + — - Pi) = Cimi . 



( 2 . 2 ) 



Summing over the entire inventory, the average on hand invest- 
ment would be 



n Q . n 

i', + -7 - >'i> ' .1 • 

1=1 1=1 



(2.3) 



Hadley and Whitin [1] also show;' the average number of 
orders (d^) per unit time for the ith item in steady state is 



X. 

= d. , 

Q. 1 



(2.4) 



where 

K cx> 

Xi = /q X f^(x)dx and fj^ (x) 

is the density of demand per unit time. Then the average num- 
ber of orders per unit time for the entire inventory would be 
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(2.5) 



n A . n 

I = I a • 

ill iSi ■- 

Therefore the constraints on average investment and on 
the number of orders can be written as 



n Q. 

I c.(r. + ^ - p ) < K 
1=1 



(2.6) 



and 



n A . 

i fT 1 ‘^2 • 

X 



i^i Q 



(2.7) 



2.2.2 The Objective Function 

Several possible objectives can be proposed. At this 
point, however, we will select as the objective function the 
minimization of expected units short per unit time. The solu 
tion technique would be similar for any choice of objective 
function. 

From Hadley and Whitin [1] we see the expected shortages 
per procurement cycle for the ith item can be expressed as 

(x-r ) f . (x,T . ) dx . (2.8) 

, r . 1 1 ' 1 

1 

The expected shortages per unit time for the ith them become 



A . 



X /“ (x-r^) f^ (x,T^) dx 



(2.9) 



1 1 

The objective can now be stated as the minimization of 
the total number of shortages per unit time for the entire in 
ventory. The objective function becomes 



Z = 



n A . 
r X fo® 



Q- -^r. (x-r . ) f • (x , T . ) dx 

1=1 '^i 1 111 
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The basic continuous review formulation can now be stated 



as 

n X ^ 

min Z = I /“ (x-r^) (x,T . ) dx , 

1=1 ^i i 



subject to 



n 

I c^(r^ + 
i=l ^ ^ 



Q 



- u^) 



1 ^1 



n X 
1=1 1 



Q. >0, and r. unrestricted. 
1 — 1 



2.3 PERIODIC REVIEW FORMULATION 

If we are again given an n item inventory and given a 
time period T between reviews and a limitation on average in- 
vestment, the periodic review model can be formulated which 
minimizes expected units short per unit time. 

2,3.1 Constraint Formulation 

From Hadley and Whit in [1] we see the average inventory 

for the ith item can be expressed as 

XT 

Ri + Ui + -^ = bf , (2.10) 

where 

R^ = the "order-up-to" level of the ith item, 

/ OO 

xh . (x , T . +T) dx where h (x,x.+T) is the demand 

distribution over a lead time plus a review period, 
Xj^ = mean demand per unit time for the ith item, and 

T = the length of time between reviews. 
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If we assume T is fixed and the same for all items, the 
(i=l , 2 , • • • ,n) are the only decision variables. 

The average investment for the entire inventory then 
becomes 



n 

I c 

i=l 



n 



X .T 

. (R. - y . = I 

11 1 2 ' . ^ 



i=l 



c. b > 
1 1 



( 2 . 11 ) 



2.3.2 The Objective Function 

From Hadley and Whitin [1] the average number of back 
orders incurred per unit time is 



1 ^ 00 

2- = m It. (x-R. ) h . (x ,t . +T) dx . 
1 T-’R. i i 'i 

1 



( 2 . 12 ) 



Therefore, the periodic review problem given a fixed re- 
view period may be stated as 

z = y - /' 

T •'] 



subject to: 



. - R. (x-R. ) h . (x ,T .+T) dx 

1=1 1 1 1 ' 1 



n X .T 

y c. (R. - y . ^) < Kt . 

11 '^1 2 —1 

1=1 
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CHAPTER 3 



SOLUTIONS 

3.1 SOLUTION TECHNIQUE 

The technique we will use to solve the problems formula- 
ted in Chapter 2 is known at the Kuhn - Tucker Theorem. This 
theorem is often used to solve constrained minimization prob- 
lems where the constraints are inequalities and/or equalities. 

First, the Kuhn - Tucker theorem [2] will be stated as 
it applies to these problems. Consider the problem; 

minimize f(X,Y) 

subject to; 

(X,Y) _< 0, i=l, • • • ,n , 

g^(X,Y) = 0, i=n+l,***,m , 

V, >0, and x unrestricted. 

■^1 — 1 

Theorem 1 . If the constraint qualifications are satis- 
fied for the minimization problem, then for (X*^,Y*^) to 
be an optimal solution it is necessary that (x'^,Y*^) and 
some ^ satisfy conditions (1) and (2) for 

m 

F(X,Y,1I) = f(X,Y) - I n.g-(X,Y) 

1=1 ^ 

and 

_< 0 , 1 = 1 , • • *n, 

ri^ unrestricted for i=n+l , ♦ • • ,m, 

(1) a. F^(X°,Y°,H°) = 0 

b. F (X Y°,1I°) > 0 and Y°F (X°,Y°,H°) = 0 

y / y 
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(2) a. _< 0 for i=l,**‘n and £ 0 



• • • 




b. (X°,Y*^,1T) = 0 for i=n+l,***m unrestricted. 



Notice the theorem does not guarantee a solution. In- 
stead, it says if we can find a solution satisfying condi- 



tions for a minimum. 

3.2 SOLUTION TO A SIMPLIFIED CONTINUOUS REVIEW MODEL 

Before solving the general continuous review model stated 
in section 2.2, we shall consider a less complicated model, 
a form of which is used in Navy inventory management today. 
This program was referred to in the introduction as the VOSL 
program. 

Suppose the order quantities are fixed by some other cri- 
terion. Specifically we assiame 

1. Order quantities are determined from the economic 
order cost equation 



where A is the order cost and I is the holding cost, 
which we assume are the same for all items, 
n X . 

2- I 0^ = K 

i=l “i 

As stated in section 2.1, the general continuous review 
problem is 



tions (1) and (2) then (X*^,Y*^) satisfy the necessary condi- 




1 



min Z 




(3.1) 
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subject to 



? Q,” 

y c,(r. + — :r--y.) <K, 
i=l ^ ^ ^ ^ ”■ 1 

n X . 

. Q. — 2 

1=1 ^1 



(3.2) 



(3.3) 



Immediately we see that assumption 1 fixes the value of 
and therefore will eliminate the second constraint in the 
general problem. While we do not actually know the ratio of 
A to I , this ratio is implied by the constraint on the number 
of orders. Assumption (1) implies that 



Q. = i 
1 K c , 
1 



(3.4) 



where 



1 /~2a~ 

X . 

However, since the quantity — is required in equation (3.3) 

Ui 

we write (3.4) as 



X . 

^ = K/c.A. . 

Substituting (3.5) in equation (3.3) yields 



(3.5) 



n 



y K/c. X. = K- 
. ^ , 11 2 
1=1 



or 



(3.6) 



K. 



K = 



n 



y /Fir 

i=l ^ ^ 

The determination of K then fixes the order quantities 
from equation (3.4) and eliminates one set of decision 
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variables from the problem. If we substitute (3.4) into equa- 
tion (3.2) , we obtain 



n 




• 

f— 1 

v| 




y c.(r. 
iil ^ ^ 


+ — i- 

2K 


(3.7) 


We can reduce (3. 


7) to 


the form 




n 


1 5 


n 




y c.r. < Kt - 

iii 1 1 - 1 


2K . ^ 
1 = 


/A .c. + y c.p. = K' . 

1 11 lii i"i 1 


(3.8) 



Our simplified problem can now be stated as 

n ^ 

min Z = y k/xTcT f°° (x-r . ) f . (x.t . ) dx 

ii-'r. i i 'i 

1=1 1 

subject to: 



n 

J c.r. < K' 

iii 

3.2.1 Necessary Conditions 

To solve the simplified problem we apply theorem 1 from 
section 3.1. From condition (1) we have 



-K/c. A . /” f.(x,T.)dX - rinC. = 0 . 

ii-'r. I 'l li 



(3.9) 



From condition (2) we have 

n 



y c.r. - K' < 0 
. ^ , 11 1 — 
1=1 



(3.10) 



and 



n 



Hi ( I c^r. - Kp = 0 . 

1 i=l 1 1 1 



(3.11) 



Observe that equation (3.9) can be solved for yield- 



ing 



^1 = 



^ / A.c. 

-K/ 1 1 



F, (rj 



c . 1 1 

1 



(3.12) 
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where 



F^^(r^) is defined as /"* f.(x,T.)dx. 

X 

The right-hand side of equation (3.12) is always less 
than zero unless F^(r^) = 0. The case where F^(r^) = 0 will 
occur if r^ exceeds the largest lead time demand. If we 
assume this is not the case, (3.10) and (3.11) reduce to 

n 

y c. r . = K' . (3.13) 



Thus the conditions for solution to the problem are: 



^1 = 



-k/i^ c, , 



(3.14) 



and 



I 



i=l 



c , r . 



1 1 






(3.15) 



It is possible to obtain a closed form solution which 

satisfies equations (3.14) and (3.15). However, generally 

we cannot solve these equations in closed form unless 

/“ (x-r^) f^ (x,T^)dx 

1 

exists in closed form and is a relatively simple function of 



Instead let us look for a general method of solving equa- 
tions (3.14) and (3.15). If we consider equation (3.14), 
which is 



= -k/ — F*5^(r ) 



C 11 
1 



(3.14! 



we notice the right hand side has an upper bound of 0 since 



F (r ) > 0 
11 — 
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and a lower bound of 



-K 




for each i since 



F^(r.) < 1 . 

1 1 — 



This also implies 




^ Kmin( 



Suppose 6 = Kmin ( ^ j for all i. Then a suggested solu- 




i 



tion procedure would be to begin at = ■j, solve equation 



(3.14) for the vector r and compute the value of the con- 
straint using equation (3.16), which is 



Compute the value of the constraint using equation (3.16). 

If the decrease (or increase) of has not caused the sense 
of the inequality to switch, decrease (or increase) by the 



solving for the vector r at each value of and computing 
the value of H from (3.16) until the sense of the ineauality 
switches again. Continue until H = or until H is within 
some acceptable region of Kj^. This method will converge to 
the optimal solution rapidly. 



n 




6 6 

If H > K| decrease t ]^ by j. If H < K^, increase by j. 



same amount j. If the sense of the inequality has changed 
then reduce the increment to -g and decrease (or increase) 
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This approach is feasible on a high speed computer and 
takes little time for a large inventory. The limitation of 



this method is solving for r^ from the equation 





However, approximations are available for some distributions 
which cannot be solved in closed form. 

3.2.2 Sufficient Conditions 

From the Kuhn - Tucker paper [2] , if we have a convex 
objective function and a convex constraint region, the neces- 
sary conditions are also sufficient. The condition for apply- 
ing the Kuhn - Tucker theorem is that the constraint region 
be convex. Since the constraint under consideration is lin- 
ear in r, the region is convex. To show Z(r^) is convex let 
us consider the equation of the expected units short per unit 
time , 



/ oo 

(x-r- ) f.(x,T . ) dx 
r . 1 i 1 



Now if 




^ 0 for all r^, then is convex. Taking partial 



derivatives , we find 




and 




(3.17) 
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Equation (3.17) will always be greater than or equal to 
zero. Under these conditions is convex. Thus it follows 
that Z is convex since it is the sum of convex functions. 

3.3 PERIODIC REVIEW, FIXED PERIOD SOLUTION 

In this section we shall consider the following problem 
as stated in section 2.3: 



min Z 



n 1 ^ 

I ^ /r (x-R^)h. ( x,t.+T) dx 
i=l i 



subject to: 



n 

I 

i=l 



c . (R . - U • - 
11 1 



A .T 
1 , 



^^1 



3.3.1 Necessary Conditions 

Applying Theorem 1 we obtain the following equations 
which will yield the necessary conditions for optimality: 

-H?(R^) = nc^T (3.18) 



and 



n 



i=l 



c . (R. - y . - 

11 ^1 



X ,T 

= K^. 



(3.19) 



The solution of equations (3.18) and (3.19) can be found 
using the method described in the previous section. If the 
solution vector R° satisfies the conditions outlined in sec- 
tion 3.2.2, then equations (3.18) and (3.19) are both neces- 
sary and sufficient. 

3.4 GENERAL CONTINUOUS REVIEW 

Let us now consider the continuous review problem stated 
in section 2.2,2 which was: 

n A. ^ 

minimize Z = T / (x-r )f (x,T.)dx , 

^ 1 O •' r o 1 1 1 

1=1 ~i 1 
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svabject to; 



n o . 

-T - «i> i ■'l ' 



n X , 

i£l Qi 1 ^2 ' 



Q. >0, and r- unrestricted . 

1 — 1 

3.4.1 Necessary Conditions 

As a result of condition (1) in section 3.1, we have the 
following equations for all i 



^ F=(r.) - n^o. = 0 
^1 



(3.20) 



and 



X. ^^2^1 

- — o /„ (x-r, )f . (x,T. )dx - — ^ + s- > 0 . (3.21) 

Q .2 111 2 q ^2 - 

Also, because of condition (1), we know (X° ,Y° , H) = 0. 

If we modify equation (3.21) by multiplying through by 0^, we 
obtain 



Hi c. 0. hoX . 

1 11 . 2 1 



O. 



- ^ /"^ (x-r.)f. (x,T^)dx - 

As a result of condition (2) we have 
n Q. 

.1 + -5 - “i> - 1 0 

1=1 



= 0 . (3.22) 



and 



n Q. 

=i<‘'i + 

1=1 



= 0 . 



(3.23) 



(3.24) 



Equation (3.20) can be solved for which yields 



'’l = - sV ''i<“'i> • 

11 



(3.25) 
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Since > 0, > 0 and 0, is always less than 

zero unless exceeds the greatest possible lead time demand. 
Assuming this is not the case, equations (3.23) and (3.24) 
reduce to 



In addition, from condition (2) we have 

n X . 

I 



(3.26) 



0 - T C . 

1 = 1 1 



and 



n X c 

Z 



. - c . 

1=1 1 



= 0 



Then considering equation (3.22) and substituting for equa- 
tion (3.25), we have 



X. 



X 



“ rT^ /°° (x-r.)f. (x,T.)dx + — ^ F?(r.) + 

Q.-'r. i i 'i 2 11 



^2^i 

Q,‘ 



= 0 . 



Solving for ri 2 yields 

Q. 



rio = - — ^ F?(r,) + /” (x-r.)f. (x,x.)dx 

2 2 11 ■'r. i i 'i 



or 



Q. 

rio = - (r. + -i)F?(R.) + /” xf.(x,T.)dx 

2 1 2 i i ■’r. I 'l 



(3.27) 



The right-hand side of equation (3.27) is not always neg- 
ative. However, theorem 1 requires n • <0. This suggests two 



1 — 

possible cases for consideration. 

Case I . ri 2 < 0, < 0 

As a result of condition (2) , the necessary conditions in 
this care are : 

- ^ F^(r,) - niCi = 0 , (3.20) 

~1 
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= 0 , (3.22) 



- ^l^jQj ^ ^2^i 



“ ^ ^r. - 



'1 1 



0 . 



n 



Q 



C.(r. + = 0 , 



(3-26) 



and 



n X . 

.1 0 ^ - ^2 = » • 
1 = 1 ~1 



(3.28) 



Case II. 



Tl2 = 0, rii < 0 



The necessary conditions in this case are 
X 



■ Q~ ■ "l^i = ° ' 



(3.20) 



- ^ (x-r^)f^(x,T^)dx " = 0 , 



1 1 



(3.27) 



and 



n Q . 

I c (r. + ^ - y ) = K . (3.25) 

i^l 1 1 z 1 1 

Since the equations above cannot be solved for r° and 
Q° in general, some iterative process or search technique 
seems desirable. One such technique was used in solving an 
example problem stated in section 5.2. 

Beginning with the solution to the simplified continuous 
review problem, which would be a feasible one, compute the 
value of T]^ for all i. Find 

h = min [ri 2 (i) ] 

and 

f = max [ri 2 (i) ] 



27 



It is now desirable to change the feasible solution ob- 
tained previously so that h increases and f decreases until 
h = f. When this occurs, we have satisfied the necessary con- 
ditions for optimality. If ri 2 (j) ~ ^ Ti 2 (k) = f, we re- 

duce Q . by 1 and increase Q, using equation (3.28) such that 
D ^ 

(3.28) is still satisfied. Then we solve another simplified 

continuous review problem using the new values of Q. and 0, . 

3 

We continue until h = f. 

The algorithm described above is an elementary approach 
to the problem. More sophisticated approaches are available. 
Among these is a search technique proposed by Fiacco and 
McCormick [ 3] . This search technique does not require a convex 
objective function or a convex constraint region. Of course, 
any solution is strictly a local minimxam in this case. The 
Fiacco - McCormick search is limited by the requirement to 
compute the inverse of a (2n + 2) x (2n + 2) matrix. This 
limits the size of the inventory for which this technique is 
useful . 

3.4.2 Sufficient Conditions 

If the objective function was convex we could apply the 
principle used in the preceding models to determine suffi- 
ciency conditions. However, Veinott (4] states that the ob- 
jective function, equation (2.8), in general is not convex. 
Veinott also discusses the conditions under which it is con- 
vex The following sufficiency conditions do not depend 
upon convexity of the objective function. 
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Consider the problem 



min f (X,Y) 

subject to: 

(X,Y) £ 0 i=l, • • • ,m 

The set of m equations g^(X,Y) £ 0 can be written as a 
vector equation 

G(X,Y) £ 0 . 

Partition G(X,Y) £ 0 such that 

G(X°,Y°) = [G^^^X°,Y°) ,G^^^ (X°,Y°) ] 

where 

g(1) (x°,y°) = 0 



and 



g(2) (x°,y°) < 0 , 



and either may be empty. 

Partition G^^^(X°,Y°) as 



( 1 ) 



** 



(X^,Y") = tG (X^,Y"),G (X^,Y^)] 



** 



where the ri . associated with each element of G (X ,Y ) is 
1 

^ o o ^ 

zero and the associated with each element of G (X ,Y ) is 
different from zero. 

Define D as the 2n x m matrix 



D = 



3g^ 3g. 

3x. ' 9y. 

1 -^1 



Partition D° and as 



D° = [D ,D ] , 

★ 



= 



H 

0 



If we define F (x,y,1I) as 



F(x,y,H) = f(x,y) - G (x,y)H 
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and define [E°] as 

r 



[E°] = 



King [5] states and proves the following theorem re- 
lating to sufficient conditions. 

Theorem 2 . In order that f (X,Y) have a local minimum 
at (X,Y) , it is sufficient that (in addition to Theorem 1) 

(X-X°,Y-Y°)'^[E°] (X-X°,y-Y®) > 0 
r 

for all X and Y satisfying 

(X-X°,Y-Y°)D* = 0 . 

It appears that Theorem 2 provides sufficient conditions for 
Case I and Case II. The only difference being, in Case I, 
that 

G (X°,Y*^) includes g^(X°,Y°) and . 

In Case II, G (X°,Y°) only includes gj^(X°,Y°) 

While Theorem 2 is interesting from a theoretical stand- 
point, it does not provide a direct comoutational verifica- 
tion of sufficiency for any particular solution. However, 
in Section 5.2 such a procedure is indicated. 




9y 3x ' 9y . 9y , 
^3 1 ^3 1 
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CHAPTER 4 



EXTENSIONS OF THE GENERAL PROBLEM 
4.1 A PROCUREMENT BUDGET CONSTRAINT, CONTINUOUS REVIEW 

Suppose instead of a limit on the average investment, 
we are given a limit on the amount of money we may obligate 
for procurement in a period of time. This constraint can 
be stated as 

n 

y c.n,Q. < B (4.1) 

iii ^ " " - 

where n^ is the nxamber of orders for the ith item per unit 
time . 

If we consider a continuous review system with back- 
orders allowed, the value of n. for the ith item is the non- 
negative integer which satisfies 

r. <a. -x.+n-Q. <r+0. (4.2) 

111 11 — 1 •'! 

where a^ = the asset position of the ith item at the begin- 
ning of the period, and 



= the demand random variable per unit time. 

Assuming the inventory position at the end of the period 



is r . + 
1 

valued. 



— j and relaxing the requirement that n be integer 



then equation (4.2) becomes 



a = - X + n Q . = r - + 
1 1 11 1 



O. 

1 



(4.3) 



Solving for n^ , 



n O 
11 




X. 

1 



and 



n = 
1 



r + 
1 



Qi 



2 



+ X 



0 



- a . 



1 



(4.4) 
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By substituting (4.4) into (4.1) we have 



n Q . 

y c. (r. + — ^ + X. - a.) < B . 
i=l 11 2 1 1 — 



(4.5) 



By taking the expected value of (4.5) , the random var- 



iable x^ becomes the parameter or expected annual demand. 
This results in the following equation: 

0_. n n 



.i, ^ * .1, =ih - I, ‘ ® • 

1=1 1=1 1=1 



(4.6) 



This constraint can be rewritten as 



1=1 1=1 1=1 



or 



n O , 

I c^(r^ + ^) < B' . (4.8) 

i=l 

Notice the constraint (4.8) is of the same form as equa- 
tion (2.6), the first constraint associated with the general 
continuous review model. From this analysis we conclude a 
constraint on obligation authority presents a problem similar 
to the general continuous review problem. 

An example of a system operating under a procurement 
budget is the Navy UICP program discussed in the introduction. 
If we desire to minimize shortages, the model is 

n X . 



min Z 




(x-r^) f^ (x,T^)dx 
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subject to: 




n 



I 




4.2 A WEIGHTED SHORTAGES FORMULATION 

The basic formulation of the problem may be extended to 
several situations. If we let d^ represent the cost of a 
shortage of the ith item, the objective can be stated as 



The addition of the weighting factor d. does not increase 



the difficulty of solving the general problem. 

4.3 MINIMIZE TIME WEIGHTED SHORTAGES 

Recently in the Naval Supply Systems Command, the num- 
ber one objective of Naval logistics management was revised to 
be the minimization of time weighted shortages. In other 
words, they desire to minimize the amount of time a customer 
must wait to receive his material. If we assume the distri- 
bution of lead time demand for the ith item is normal (]i^,cr^) 
and a continuous review system, the objective function can 
be stated, Hadley and Whitin [1] , as 



min Z 




n 

min I B(Q,r) 
i=l 



where 



B(Q,r) = §t3(r) - 3(r+Q)] 



and 



3(r) = i[a^+(r-y)^]F^(r) - |(r-y)f(r) . 
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We have formulated the problem since this problem is one 
of primary interest to the Navy today. However, no attempt 
was made to provide a solution procedure or to determine 
necessary and sufficient conditions for optimality. 
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CHAPTER 5 



SOME EXAMPLE PROBLEM SOLUTIONS 
5.1 AN EXAMPLE OF THE SIMPLIFIED CONTINUOUS REVIEW MODEL 

Let us consider an inventory of three items. The solu- 
tion technique for n items would be similar. We assume the 
distribution of lead time demand is normal with mean u. and 

variance a for the ith item. These items have the follow- 
1 

ing characteristics. 



Parame ter 


Item 1 


I tern 2 


Item 3 


^i 


1000 


1500 


2000 


c , 
1 


1 


10 


20 


^1 


100 


200 


300 


a - 
1 


100 


100 


200 


II 

CO 


000 


^2 = 15 


. 



The problem is then 

3 

min Z = I K/c . A j°° (x-r , ) f^ (x,t^) dx 
1=1 ^ ^i 

subject to: 

3 Q. 3 

I c. (r^ + -y) 1 sooo + I . 

i=l ^ ^ ^ i=l 

In section 3.2, equation (3,6) determines the value of 
the constant K which, in turn, fixes the order quantities. 
Initially, then, the value of K is determined by 



35 



K, 



K = 



n 



= .04244 



y /5TX7 

i£i ^ ^ 



utilizing equation (3.4), 



'A. 

1 



, the order quantities 



can be determined as 



Q, = 746 , 



Q 2 = 289 , and 
Q 3 = 236 . 

Since the have been determined, the problem now can 



be written as 



min 



Z= y .04244/c~rr I” (x-r,)f. (x,x.)dx 

1 1 ■’ r^ 1 1 ' 1 



subject to: 

3 3 3 0 . c . 

I c^r^ < 8000 + I c.y. - I — ^ = 11,922 . 

i=l i=l i=l 

From equation (3.12) we see the solution to the above 

problem is that vector r° such that 

nc. 



F^(r ) = -— ^ 

^ k/ITc^ 



, for all i . 



(3. 12) 



1 1 



For each value of n there will be some vector r for 
which equation (3.12) is satisfied. However, since the ob- 
jective function is convex, there exists only one vector r 
such that equation (3.12) and the constraint (3.8) are satis- 
fied simultaneously. As we decrease n from zero, we see that 
at n =-. 102 , 

F^(rj^) = .0760 , 



36 



.1965 , and 




F3(r3) = .2410 



which implies 



= 243.30 



= 285.40 



and 



= 440.80 . 



Checking the constraint, we find 



3 

y c.r. = 11923.30 
i=l ^ ^ 



which is within two dollars of the required average invest- 
ment limit. 

The order quantities and reorder points determined yield 
an expected number of shortages per unit time for each item. 
The expression for the expected number of shortages per item 
per unit time is 



When the distribution of lead time demand is normal, 

(5.1) can be written, Hadley and Whitin [1], as 



Using (5.2), the expected shortages per unit time for 
each item is computed as 




(5.1) 




(5.2) 



where 




37 



^ [- (143. 30) (.0760) + 100(.1428)] 



4.52 



/ 



Z 



z 



2 



1500 

289 



[- (85.40) (.1965) + 100 (.2770)] 



56.7 , and 



Z 



3 



2000 

236 



[-(140.8) (.2410) 



200(.3114) ] = 240.0 . 



The total expected shortages per unit time for the in- 
ventory under consideration is 
3 

Z = I Z. = 4.52 + 56.70 + 240.00 = 301.22 . 

dl ^ 

5.2 GENERAL CONTINUOUS REVIEW EXAMPLE 

Once again let us consider the inventory of 3 items from 

5.1. We shall again assume the distribution of lead time de- 

2 

mand is normal (y.,a. ) for the ith item. The items under 

1 1 

consideration have the following characteristics. 



Parameter 


Item 1 


Item 2 


Item 3 


X = 

1 


1000 


1500 


2000 


C 0 
1 


1 


10 


20 


^i 


100 


200 


300 


a 0 

1 


100 


100 


200 



= $8,000 K 2 = 15 . 



Reviewing section (3,4) we note that the solution (r ,0 ) 
must satisfy 

- ^F=(r,) = niCi , (3.18) 



■ 0^ li. (x-r^)f^(x,T )dx - 
1 1 



’^2^i 



= 0 



(3.20) 



38 



(3.23) 



and 



i=l 



Oi(r. 



0 . 

1 

2 “ 



K, 



3 X . 

I = K 
i=l ^ 

Equation (3.18) is now rewritten as 



(3.24) 



rij_ = - — F?(r.) , for all i . (5.3) 

i^i ^ 

Substituting (5.3) for rij_ in (3.20), we get 

r \-2 = - J F . (r.) + / (x-r.)f. (x,x.)dx . (5.4) 

^ XTfl XX X 



In other words , the necessary conditions for solution 
are present when 

''l = - 

11 J~j 

for all i and j , and 

02 = /~ xf^(x,T^)dx - (j + r^)F^(r^) 

= / xf (x,T,)dx - (i + r.)F?(r,) , 

Jrj j ' j' '2 3 3 1 ' 



for all i and j. In addition, both constraint equations 
(3.26) and (3.28) must be satisfied. 

Using the search routine described in section (3.4.1), 
this example was solved with the following results. 



Variable 


1 


2 


3 


r 


269.65 


307.11 


409.79 


Q 


483.08 


229.12 


313.31 


'^l 


-.0930 


-.0930 


-.0930 


^2 


-9.013 


-9.012 


-9.012 
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The value of the objective function (Z) was 285.20. Since 
< 0 and ri 2 < Of Case I applied and Case II was not appli- 
cable. If ri 2 had been non-negative, we would have used the 
conditions of Case II for solution. 

Next we must determine whether the solution obtained 
above satisfies the sufficiency conditions stated in theorem 
2 in section 3.4.2 for a relative min. That is we must 
determine if 



h^Eph > 0 , 



for 

hD* = 0 . 



* 

E , h and D are defined as in section 3.4.2. 

r 

The theorem stated below is indicated by Hadley [6] and 
proved in Hancock [7]. 

Consider the matrix A 



A = 



- XI 

* 

D 




Theorem 3. If the roots of j iA| | are all positive then 
f(X,Y) takes on a strong relative minimum at (X°,Y°) . 

The procedure for applying theorem 3 is simple enough. 
Construct the matrix A. Take the determinant of A and set 
it equal to zero. The result will be an (2n x 2) degree 
polynomial in X. Using Decartes rule of signs, the number 
of positive roots of | lAj j can be determined. For example, 
suppose ||a|| =A+2X+1=0. Here we have no sign 
changes. According to Decartes rule of signs there are no 
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2X + 1 



0 



2 

positive roots. However in the equation X 
there are 2 sign changes which indicates 2 positive roots. 
To show sufficiency of the solution states previously 



we form the i 


matrix A 


evaluated 


at (Q° 


,r°) which 


is 








.196-X 


0 


0 


000193 


0 


0 


1 


0 




0 


1.475-X 


0 


0 


.00405 


0 


10 


0 




0 


0 


2.19-X 


0 


0 


. 006 


20 


0 


A = 


.000193 


0 


0 


000126 


0 


0 


.5 


.00428 




0 


.00405 


0 


0 


-.000696 


0 


5 


.0286 




0 


0 


.006 


0 


0 


0018 


10 


.0205 




1 


10 


20 


.5 


5 


10 


0 


0 




. 0 


0 


0 


00428 


.0286 


0205 


0 


0 



Taking the determinant of A and setting it equal to zero 
we find all of the roots are positive indicating the solution 
(Q°r°) is a strong relative minimum. 

5.3 SIMPLIFIED CONTINUOUS REVIEW WITH A PROCUREMENT BUDGET 
CONSTRAINT 

In this section we shall consider a three item inven- 
tory and again assume the distribution of lead time demand 
for the ith item is normal (y^/O^). The items under consid- 
eration have the following characteristics 



Parameter 


1 


2 


3 


X. 


1000 


1500 


2000 


1 


c . 


1 


10 


20 


1 


^i 


100 


200 


300 


Oi 


100 


100 


200 


a . 


100 


200 


300 


1 


B = $64 


,000 


K2 = 


15 
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The general problem is 



3 X 



min 



Z = I 7 ^ r (x-r. ) f . (x,T. )dx 
• I Q- r - 1 1 ' 1 



i=l -"i 



siabject to: 



3 Q. 

y c.(r. + -i+X. - a.) <B 
^,11 2 1 1 — 



i=l 

3 X. 

y -^ < K„ . 

. O, — 2 

1=1 “1 

Since X,, a., and c. are known constants, the first con- 
11 1 

straint can be rewritten as 



3 Q. 3 3 

y c. (r. + — i) < 64,000 - y X.c. + J a.c., 

1 1 2— ' 11 ii' 



or 



i=l 



i=l 



i=l 



3 Q. 

I c^(r^ + -j) 1 64,000 - 56,000 + 8100 = 16,100. 



i=l 

If we also assume 



(5.5) 



O 






^1 K C 



as was done in section 3.2 and section 5.1, equation (5.5) 
reduces to 

n 

1=1 1=1 



n Q . 

y c.r. < 16,100 - y c. — ^ , 

.^,11— ' i2' 



(5.6) 



or 



n 

y c.r. < 11922 . 
i=l ^ ^ - 

The problem can now be stated as 

3 X. 



min 



^=1,0^ C. (x-r,i)r (x,T.)dx 
1=1 1 
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subject to: 



n 

y c.r. < 11,922 . 

iil ^ ^ - 

Observe, however, that this problem was solved in section 

5.1. 
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CHAPTER 6 



SUMMARY AND CONCLUSIONS 

It is apparent that the models proposed do not represent 
the ultimate answer in multi-item inventory theory. The sim- 
plified continuous review problem stated in section 3.2 rep- 
resents a model of the Navy's VOSL program. From the example 
problems we see the general problem provides a solution which 
is better than the solution of the simplified continuous re- 
view problem. The major advantage of the simplified contin- 
uous review model is its computational ease. 

While the solution to the general problem as stated is 
theoretically correct, an efficient algorithm for rapid loca- 
tion of stationary points and subsequent verification of 
sufficiency conditions is at present not available. There 
is certainly room for future research on this topic. 

The assumption that the expected number of orders which 
can be processed is a well known constant is perhaps optimis- 
tic. However, if we look at the general problem as a two 
criterion problem, we can generate the following efficiency 
curve. One method of generating this curve would be to 

* 

Z * 

‘ / 

I / 

/f 
/ , 

/ 1 
i 




Figure I K 2 
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solve the general problem for several values of the parameter 

K 2 (the expected number of orders) . If these values of K 2 

versus the expected shortages per unit time are plotted, the 

curve represented by Figure I is generated. 

It is apparent that for every organization of N people, 

there is some maximum number of orders which can be processed 

* 

per unit time. This point is represented by K 2 . However, 
we contend that each individual manager must examine the 
alternatives represented by the efficiency curve in Figure 
I and select that point at which he desires to operate. 

Another approach to the problem would be to formulate 
a vector minimization problem. For example, 

min [Z,K 2 ] ; 



subject to: 



3 Q. 

y c . (r , + —TT - u . ) < K. . 

11 2 ^1 — 1 



Hadley [6] indicates a procedure for solving problems of 
this type. Such an approach to the problem proposed repre- 
sents a fertile area for future work. 

From the example problems given, it appears the Navy's 
method of fixing Q and solving for r^ (as we did in the sim- 
plified formulation of section 3.2) does not result in an 
optimal solution to the problem of minimizing shortages per 
unit time. The use of the unknown parameters, holding cost, 
shortage cost, and order cost presents a possible source of 
error. 
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We contend the constraints on average investment and 
the expected niJinber of orders actually imply the values of 
the unknown parameters. Since the values of the constraints 
are more easily determined than the order cost and holding 
cost, the models proposed seem much more appropriate than 
the traditional variable cost minimization models. 

While an efficient algorithm for solution of the gen- 
eral problem has not been presented, the advent of high 
speed computers has opened the field of iterative solution 
procedures. It should be only a matter of time until a pro- 
cedure is available which can be reasonably applied to a 
large inventory. However, the techniques discussed in sec- 
tion 3.4.1 are feasible only for small inventories or sub- 
sets of the larger inventory. For instance the problem 
could be solved for the entire inventory using the simpli- 
fied continuous review model. As stated previously, the 
simplified continuous review algorithm is computationally 
feasible for large inventories. The inventory manager could 
then select subsets of items whose decision variables in- 
tuitively appear to be unreasonable. We can then formulate 
and solve the general model for the subset of items using 
the budget allocated, to those items by the simplified con- 
tinuous review model as The number of orders allocated 

to the subset of items would become the constant K 2 . Selec- 
tion of a subset of the n item inventory is necessary since 
the general solution algorithms mentioned are feasible only 
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if the number of items is small. We can guarantee that the 
value of the objective function will at worst be the same 
as the simplified continuous review model and in all likeli- 
hood, it will decrease. 

While the multi-item problem has been solved when funds 
are unlimited, the assumption must be made that there are 
no interactions among items (i.e., the problem degenerates 
into N single item problems). For instance, it must be 
assumed that enough materials handling equipment is avail- 
able to handle all material, enough warehouse space is avail- 
able and numerous other possible interactions do not exist. 
Each of these interactions, including a limitation on funds, 
represents a constraint on some resource within the system. 

It seems logical, then, that the next step in the formula- 
tion of multi-item inventory problems should be of the form 
presented in this study. 
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